Gaussian approximations for the exchange-energy functional of current-carrying 
states: Applications to two-dimensional systems 
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Electronic structure calculations are routinely carried out within the framework of density- 
functional theory, often with great success. For electrons in reduced dimensions, however, there 
is still a need for better approximations to the exchange-correlation energy functional. Further- 
more, the need for properly describing current-carrying states represents an additional challenge for 
the development of approximate functionals. In order to make progress along these directions, we 
show that simple and efficient expressions for the exchange energy can be obtained by considering the 
short-range behavior of the one-body spin-density matrix. Applications to several two-dimensional 
systems confirm the excellent performance of the derived approximations, and verify the gauge- 
invariance requirement to be of great importance for dealing with current-carrying states. 
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I. INTRODUCTION 

The success of density-functional theory [l[ (DFT) cru- 
cially depends on the availability of good approxima- 
tions for the exchange-correlation (xc) energy functional. 
Great progress has been achieved beyond the commonly 
used local (spin) density approximation [L(S)DA] by 
means of, e.g., generalized-gradient approximations, or- 
bital functionals, and hybrid functionals [2j. However, 
most of the approximations developed so far have focused 
on three-dimensional (3D) systems. 

Most density functionals developed for 3D perform 
poorly when applied to two-dimensional (2D) systems [3j, 
0i IE @| • However, at present 2D structures constitute 
a large pool of applications in semiconductor nanotech- 
nology, e.g., quantum Hall systems, spintronic devices, 
and quantum dots (QDs) Q. Within DFT these sys- 
tems are often treated using the 2D form of the LSDA 
employing the analytic expression of exchange energy 
of the 2D electron gas [|[ (2DEG) and the correspond- 
ing correlation energies parametrized via quantum Monte 
Carlo calculations first by Tanatar and Ceperley @ and 
later, for the complete ran ge o f collinear spin polariza- 
tion, by Attaccalite et al [10|. Despite the relatively 
good performance of the LSDA in terms of total ener- 
gies, there is still a clear need for accurate density func- 
tionals in 2D. Only recently, xc functionals tailored for 
2D have been pushed forward [H El El El EE El]- 
Here we make another natural step along these directions 
by studying Gaussian approximations (GAs) for the ex- 
change energy functional. We extend previous studies on 
GAs [HI EE El, HE SI E2 to (i) 2D systems and to 



(ii) current-carrying states for both 2D and 3D systems. 
Numerical tests of the resulting exchange energies show 
significant improvement over the LSDA. 



II. GENERAL FORMALISM 

Within current-spin-density functional theory [27], HH 
(CSDFT), the ground state energy E v , the spin densities 
p CT (r) and paramagnetic current density j CT (r) of a system 
of N = JV-j- + JVj_ interacting electrons are determined by 
minimizing the total-energy functional which is given in 
Hartrce atomic units as 

Ev[{pM]=T,[{pM] + E a \p] 

/rf D r t)j (rK(r) + ^ / d D r A CT (r)j CT (r) 



<T=U 



Jd D rAl(r)p a (r)+E xc [{ P<J },U t7 }] (1) 



where D = 2,3 is the dimensionality of interest, 
T s [{p a , jo-}] is the kinetic energy functional of non- 
interacting particles, ^(r) is an external (local) spin- 
dependent potential acting upon the interacting system, 
A CT (r) is an external (local) spin-dependent vector poten- 
tial (note that physical fields are spin- independent, but it 
is formally convenient to introduce the spin-dependence), 
_Eh[p] is the classical electrostatic or Hartree energy 
of the total charge density p(r) — Pj(r) + Pi(r), and 
E xc [{pa, jo-}] is the xc energy functional. This quantity 
is usually further decomposed into the exchange and cor- 
relation energies as 
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E xc [{p„,ja}] 



E x [{Pa,3a}]+E c [{p a ,j a }], (2) 
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An exact expression is known for the exchange energy 



Ex[{Pa,iv}] 



d ri d r 2 —, r 1 

ri-r 2 



(3) 

where 7 CT is the one-body spin-density matrix (1BSDM) 
of the Kohn-Sham (KS) system 



r0i,r 2 ) =^^fe l0 .(ri)^fe i0 .(r 2 ) , 



(4) 



fc=i 



and ipka{ r ) are the KS spin orbitals. Note that we are 
assuming that the KS ground state has the form of a 
single Slater determinant. The number of electrons with 
spin a is given by 



J d D n J d D r'\ la (r u r') 



(5) 



Introducing the average and relative particle coordinates, 
respectively, as 



1 , 

r = 5 ( r i + r 2) > s = r x - r 2 , 



(G) 



Eq. ([31) can be rewritten as 

E x [PT,Pl}=-^ £ j d ° r j s [D ~ 2) ds(\lA 2 ), (7) 
where 



1 I 2 








\1<x\ 


^ J 


[ dtt D 


7a 



r H — , r 

v 2 2 



(8) 



is the angular average of |7 CT | 2 . Note that (|7er| 2 ) is a 
function of the position, r, but this fact is omitted in the 
notation. Also in the remaining part of the present work, 
we often shorten the notation where the dependence on r 
is clear. The term dQu in Eq. (|HJ) indicates the angular 
differential in D dimensions, Q 2 = 27r, and Q3 = 47r. In 
the new coordinates, Eq. |(SJ) can be rewritten as 



N a =Cl 



D 



d D r / s^dsd^] 2 ) 



(9) 



It is apparent that an approximation for (Ko-I ) in Eq. 
(|8]l provides an approximation for the exchange energy in 
Eq. ([7]). Equation © allows us to check the validity of 
such an approximation in terms of the particle-number 
normalization. 



III. GAUSSIAN APPROXIMATION FOR THE 
ONE-BODY SPIN-DENSITY MATRIX 

One of the more successful strategies to develop ap- 
proximation for the exchange energy functional is based 
on the density-matrix expansion scheme. [HI, [U, [25|, [26j 



Here we reconsider the so-called short-range behavior of 
the 1BSDM. In particular, we extend the results previ- 
ously derived for zero-current states in 3D [HI, [T^j to 
current- carrying states of both 3D and 2D systems. 

In dealing with current-carrying states it is necessary 
to allow the orbitals in Eq. ^ to be complex valued. 
Expanding the 1BSDM in the inter-particle coordinate, 
s, at a given average position, r, we find 



where 



•T 



1 iV„ 

-Y 

k=l 



(10) 



(ii) 



is the spin-dependent (paramagnetic) current density, 
and 



1 V"^ Jkrr.aJkiT.P 



with 

[t<T 



a, (3 



1 dpka dp ka 



■f pka dx a OX/3 dx a dxj} 



Pa 



(12) 



(13) 



where a, (3 runs over the Cartesian coordinates in 3D 
or 2D. The term [S CT ] a g can be interpreted as kinetic 
energy-density tensor. We point out that So- is not a 
gauge-invariant quantity due to the second term in Eq. 

At this point, we carry out an exponential resumma- 
tion of Eq. (JTOJ), which yields 



Pa- exp 



b t S ct s 



•T 1 



(14) 



where s T and j T are the transposes of the vectors s and 
j, respectively. Although this approximation reproduces 
the exact short range behavior of the 1BSDM, it seems - 
at first sight - a rather oversimplified approximation for 
arbitrary s. However, it has the appealing feature that 
the corresponding Wigner transform [l| reproduces the 
spin-particle density, the paramagnetic current density, 
and the kinetic-energy-density tensor of the KS system 
exactly (see below). 

The Wigner transformation of Eq. (|14|l is given by 



7W i \ 

fa (r,p) 



^/(2tt) d det (fc B © CT ) 
■ \ T 

— (k B ®o 
Pa J 



x exp 



- 5 IP- 



P 

Pa 



(15) 



where, following Ref. [Tj|, we define the matrix ® a to 
satisfy the relation 



1 



-k B Pa®c 



(16) 
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where fcs is the Boltzmann constant. The presence of 
the Boltzmann constant is motivated by the fact that for 
the spin-unpolarized zero-current case, the Wigner trans- 
form in question provides a description for the KS state 
formally resembling a thermodynamic one (l?! [l8l [T9L I20I 
I2H I22I ] . Instead, in this work we allow the orbitals to be 
spin-unrestricted and complex valued. We observer that 
& a is apparently not gauge- invariant, thus we shall re- 
frain from interpreting it as a local temperature. Finally, 
it is immediate to verify that 



with 



A,(r)= / d"p/f(r,p) 



j CT (r)= / d u pvff (r,p) 



and 

[^L )/3 (r) 



1 f jD ( 3"A V ) 

- d p [p a -— 

fW 



(17) 



(18) 



PP 



x /<T(r,p). 



JaA r ) 

Pa{?) 

(19) 



IV. GAUSSIAN APPROXIMATIONS FOR THE 
EXCHANGE ENERGIES 

In this section, we provide few approximations for 
the exchange energy of spin-polarized current-carrying 
states. Although what follows may be applied to the 3D 
case as well, we shall restrict ourselves to 2D systems 
where the need of new approximations in the quantum- 
Hall regime, for example, is particularly large. 

Let us first observe that substitution of Eq. (fT4|) in Eq. 
([8|) would yield an expression which is not gauge invari- 
ant. Popular (3D) meta-generalized gradient approxima- 
tions for the xc energy are also not gau ge in variant [29| . 
This fact is not only formally inadequate [27l. [28| , but it is 
also a source of practical problems. For example, it has 
been observed that the lack of gauge invariance causes 
a wrong description of degenerate atomic ground states 
carrying different paramagnetic currents. Recently, Tao 
and Perdew have proposed a correction for this prob- 
lem [29T |30| . Here we show how this kind of correction 
naturally emerges in the context of the present GA. 

In order to fulfill the gauge invariance requirement, we 
substitute Eq. (fTU|) in Eq. ©, and obtain 



<w 2 > 



where 



Tr(S ff ) 



1 - 



E 

fc=i 



Tr(S c 

Pa 



1 (Vp fco 



Pko 



1 (Sa_ 

2 [.Pa 



2 Pka 



(20) 



-y 2 pa 



(21) 



(22) 



fc=i 



being the spin-dependent kinetic-energy density. Note 
that, while r a is not gauge invariant, the modified quan- 
tity 



T n =T„- 



1 fSl 



2 \p. 



(23) 



is indeed gauge invariant. Now, let us perform an ex- 
ponential resummation of Eq. ([20]) . which leads to the 
following GA 

<| 7ff | 2 )^exp(-0 (24) 



where 



PA = 



iv 2 Pff 1 fi 



Pa 8 p a 2 \p 

or, alternatively, using Eq. (|16[) 

(3A = \ Tr(e r )-(k 

L \P<y 



(25) 



(26) 



Expression (|26|) shows how f3 a relates to ® a ■ In contrast 
to @<j, p a is a gauge- invariant quantity. Hence, the ap- 
proximation for the exchange energy obtained from Eq. 
(|2"6")l is gauge invariant as well. 

We further focus on inhomogeneous 2D systems. An 
expression for the exchange energy is readily obtained by 
inserting Eq. (gj) in Eqs. © and ©: 



E x[PhPl] = - '- 



r 3/2 



E 



d 2 rpl(r)A /2 (v). (27) 



This is a functional for the exchange energy in the form of 
a current-dependent meta-generalized gradient approxi- 
mation. To stress the dependency on the paramagnetic 
current, and to remind its origin from a GA, we will in- 
dicate the functional as J-GA. 

Furthermore, we impose that our approximation repro- 
duces the correct normalization condition (|5|) . Hence, we 
modify Eq. ([24l by adding a fourth-order term [2fj| as 
follows 



Pi exp 



Pa 



(28) 



This expression leads to 



E x [p hPl ] = E / 



x pl{v)Aj 2 {v). 



(29) 



4 



The parameter A a is determined by using the following 
relation 



N n = 



<fr[l 



2A a ]pl(v)l3 a {v). 



(30) 



Eq. ([29]) together with Eq. (f30j) provide another density 
functional for the exchange energy. We refer to this func- 
tional as the current-dependent modified GA denoted by 
J-MGA. 



V. APPLICATIONS 

In order to test the performance and the degree of uni- 
versality of the J-GA and J-MGA, given in Eqs. (J57J) and 
([2"9"| . respectively, we consider below different 2D elec- 
tron systems, i.e., QDs of different shapes, as well as the 
2DEG. We also assess the importance of the approxima- 
tion to account for the paramagnetic current density in 
Eqs. PD|) and (|2"Tj) , and thus the prerequisite of the 
gauge invariance. The assessment is done by comparing 
J-GA and J-MGA to the results obtained by neglecting 
in the functionals the terms depending explicitly on 
These approximations are denoted below as 0-GA and 
0-MGA. 



A. Finite systems 

As relevant examples of finite 2D electron systems 
we consider parabolic (harmonic) 0] and rectangular 
QDs [HI, respectively. In both cases, the many- 
electron Hamiltonian is given by 



H 



N 

£ 



5 I p. 



C 



+ F ext (ri) 



N 

l<j 



1 



( 31 ) 

where N is the number of electrons and A is the external 
vector potential of the homogeneous magnetic field B = 
Bz perpendicular to the 2D (x-y) plane. 

The external confining potential is defined for a 
parabolic QD, containing here N = 2 ... 20 electrons, as 



(32) 



where loq is the confinement strength. For N — 2 we have 
u) = 1 (see below), and otherwise oj = 0.42168 corre- 
sponding to a typical confinement of 5 meV when apply- 
ing the effective mass approximation (with the effective 
mass m* = 0.067 m e and dielectric constant e = 12.4 eo) 
in the modeling of GaAs QDs Q. 

The rectangular QD containing N = 4 ... 16 electrons 
is defined by 



0, < x < aL, < y < L 
oo, elsewhere, 



TABLE I: Exchange energies for closed-shell parabolic (upper 
part) and rectangular (lower part) quantum dots calculated 
using the exact exchange (EXX) (in the KLI approximation 
for N > 2), the Gaussian approximation (GA), the modified 
Gaussian approximation (MGA), and the local spin-density 
approximation (LSDA). The relative errors with respect to 
EXX are also shown. 



N 


EXX 


J-GA 


A (%) 


J-MGA 


A(%) 


LSDA 


A (%) 


2 


-1.08 


-1.12 


-3.0 


-1.10 


-1.5 


-0.98 


9.3 


6 


-2.23 


-2.28 


-2.1 


-2.28 


-2.2 


-2.13 


4.4 


12 


-4.89 


-5.01 


-2.5 


-5.03 


-2.9 


-4.76 


2.6 


20 


-8.78 


-9.00 


-2.5 


-9.05 


-3.0 


-8.63 


1.7 


6 


-3.14 


-3.33 


-5.9 


-3.25 


-3.3 


2.99 


4.9 


12 


-8.19 


-8.46 


-3.3 


-8.42 


-2.8 


-7.99 


2.5 


16 


-12.7 


-13.3 


-4.4 


-13.1 


-3.2 


-12.3 


3.5 



(33) 



where L = y/2n is the (smaller) side length, and a — 2 
determines the side ratio of the rectangle. The size is then 
~ 90 nm x 45 nm. The electronic properties of similar 
rectangular QDs have been studied in detail in Refs. 31] 
and [321 ]. 

We point out that, in principle, the presence of A in 
the Hamiltonian brings out the need of the full-flagged 
current-spin-density functional theory [27| . However, at 
the level of the KS equations we neglect the xc vector 
potential, but still include the external one. This ap- 
proximation has been shown to be accurate for atom- 
istic systems, and QDs up to relatively high magnetic 

fields [ni, m m, hm m wa m . 

As reference we use the exact- exchange (EXX) results 
calculated in the Krieger-Li-Iafrate (KLI) approxima- 
tion [3^]. The converged KS orbitals from the KLI are 
then used as the input for our functionals. An exception 
is the two-electron parabolic QD with lvq = 1, for which 
we can apply the known analytic density [40| in the cal- 
culation of the EXX energy (which is minus half of the 
Hartree energy) and as the input in the functionals. For 
comparison, we also compute the exchange energies from 
the (2D) LSDA ||. Both the EXX and LSDA calcula- 
tions are performed using the octopus real-space code 
within SDFT [il|. 

Table Q] shows the exchange energies for a set of both 
parabolic (upper part) and rectangular (lower part) QDs. 
The relative errors with respect to EXX are also shown. 
All the cases correspond to closed-shell ground-state solu- 
tions possessing no currents, and the total spin is S = 0. 
Therefore, the J-GA (J-MGA) result is similar to 0-GA 
(0-MGA). Overall, we find a good agreement between 
the J-GA and the EXX. The J-MGA, which contains the 
fourth-order term in (|7o-| 2 ), is in most examples more 
accurate: The deviation from the EXX is < 3%. As ex- 
pected, the LSDA generally overestimates the exchange 
energies, but the deviation from the EXX reduces as a 
function of N. In contrast, J-(M)GA seems to gradu- 
ally lose its accuracy when N is increased, especially in 
parabolic QDs. This may be due to the fact that since 
the electron density profile is relatively flat in large QDs, 
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the LSDA correspondingly becomes more suitable. 

In Table |TT] we show the exchange energies for spin- 
polarized current-carrying states in parabolic (upper 
part) and rectangular (lower part) QDs. The N = 2 
QD with S = 1 corresponds to an excited state, whereas 
in the other systems the spin polarization of the ground 
state is achieved by applying the external magnetic field. 
Compared to the EXX, the accuracy of both J-GA and 
J-MGA is excellent. Moreover, both functionals are con- 
siderably more accurate than the LSDA. On the other 
hand, the approximations 0-GA and 0-MGA obtained by 
omitting the current-dependent terms in Eq. (|2H)) lead to 
completely wrong results. This demonstrates the impor- 
tance of the currents involved in the derivation of a gauge 
invariant GA. 



This expression can be used as an explicit density func- 
tional in the LSDA fashion. In fact, the only difference 
to the LSDA exchange [|[ is the prefactor, which is here 
~ 4.4% smaller than in the LSDA. Therefore, one can im- 
mediately deduce from Tables HI and HT1 that . when consid- 
ering finite QD systems of few electrons, Eq. ([33)) leads 
to good agreement with the EXX results by correcting 
the LSDA exchange energies by ~ 4.4%. 

Interestingly, Eq. (|35[) has been recently obtained in an 
alternative way by considering the short-range behavior 
of the exchange-hole function This identity of the 

expression can be seen as a rather reassuring result. 



VI. CONCLUSIONS 



B. Two-dimensional electron gas 

Finally let us consider the GA for the exchange energy 
in the case of the 2DEG. We restrict ourselves to the 
case of vanishing external vector potential, but still we 
allow the external scalar potential to be spin-dependent. 
In this situation, we can set \7 2 p a = and Vp a = 0. 
Moreover, we may set j CT = by considering the 2DEG 
at rest. Hence, the current term in Eq. (|2"5|) drops out 
from the final expression 



(I7.I 2 ) 



Pa ex P 1 



2 2\ 



(34) 



We note that correct particle numbers for each spin chan- 
nel are obtained from Eq. (|34p. and therefore A a = 0. 
Defining the 2D density parameter r s = 1/ y/np, and the 
polarization £ = (p-p — Pi)/p, the total exchange energy 
per particle becomes 



4^2? 



+ 



3/2 



3/2 



(35) 



In this work we have analyzed the one-body spin- 
density matrix for current-carrying states. We devised 
gauge-invariant approximations for the exchange en- 
ergy by considering a Gaussian approximation for the 
short-range behavior of the one-body spin-density ma- 
trix, in particular, for its angular average of the module 
square. The resulting simple approximations for the ex- 
change energy perform extremely well in a variety of two- 
dimensional systems, including a diverse set of quantum 
dots, as well as the homogeneous two-dimensional elec- 
tron gas. Moreover, our analysis confirms the relevance of 
the gauge invariance when dealing with current-carrying 
states. 
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